Stability of edge states and edge magnetism in graphene nanoribbons 
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We critically discuss the stability of edge states and edge magnetism in zigzag edge graphene 
nanoribbons (ZGNRs). We point out that magnetic edge states might not exist in real systems, and 
show that there are at least three very natural mechanisms - edge reconstruction, edge passivation, 
and edge closure - which dramatically reduce the effect of edge states in ZGNRs or even totally 
eliminate them. Even if systems with magnetic edge states could be made, the intrinsic magnetism 
would not be stable at room temperature. Charge doping and the presence of edge defects further 
destabilize the intrinsic magnetism of such systems. 
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I. INTRODUCTION 

Graphene is a single layer of graphite. Since the 
first successful isolation of mono- and multi-layers of 
graphene, 1 this material has raised tremendous atten- 
tion within the scientific community. This is mainly 
due to the fact that graphene is an almost perfect two- 
dimensional system with unique electronic properties, 
that might ultimately lead to a new generation of na- 
noelectronic devices. 2,3 

Graphene nanoribbons (GNRs) are one-dimensional 
stripes "cut" from graphene, which are a few nanometers 
in width and quasi infinite in the perpendicular direc- 
tion. Initially GNRs were primarily discussed in the the- 
oretical literature^ but meanwhile GNRs can be fabri- 
cated and studied in the laboratory £~— GNRs with zigzag 
edges (ZGNRs) are of particular interest, because they 
are forming edge states^ The latter are localized elec- 
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FIG. 1. Graphical outline of this article. Top: Three mecha- 
nisms that strongly affect the edge states in ZGNRs: (a) the 
atomic structure of the 2211 passivated edge, (b) the atomic 
structure of the zz(57) reconstructed edge, (c) edge closure. 
Bottom: The edge magnetism in (d) ideal ZGNRs is not sta- 
ble at room temperature. The low temperature magnetism is 
destabilized (or even totally disappears) whenever the system 
has (e) edge defects or (f) is charge doped. 



tronic statesthat decay exponentially towards the center 
of the ribbon. 2 The decay lengths are in the range of 
a few nanometers J£ Such edge states do not only ex- 
ist in perfect ZGNRs, but in any graphene system that 
has zigzag edge segments ii The localized nature of the 
edge states in ZGNRs gives rise to flats bands (i.e. parts 
of bands with little dispersion) in the electronic band 
structure, and thus to a pronounced peak in the elec- 
tronic density of states (DOS) right at the Fermi level 
(E-p). This peak has been measured in the local DOS 
of monoatomic zigzag step edges of graphite by scanning 
tunneling microscop y 10 i 12 i 13 In those measurements the 
peak in the DOS did not appear right at Ey, but very 
close to it. However, no such measurements for a mono- 
layer of graphene or GNRs are known up to date. 

According to numerous theoretical studies based on 
density functional theory or on the mean-field Hubbard 
model^ii^ it is believed that the peak in the DOS at E? 
gives rise to a magnetic instability, where the edge states 
become spin-polarized. In the magnetic ground state 
a band gap at E? is opened, and the atoms are ferro- 
magnetically ordered along one edge, and antiferromag- 
netically ordered between opposite edges [see Fig. EJa)]. 
This antiferromagnetic ground state is consistent with 
the Lieb theorem for a bi-particle lattice within a Hub- 
bard model type of description^ Potential applications 
of edge magnetism for spintronics are widely discussed 
(for an overview see (l4(). But up to now, the edge mag- 
netism in ZGNRs was never directly observed in local 
probe microscopy. Only an indirect observation has been 
reported quite recently^ 

In this article we critically discuss the present state of 
research on graphene edge states and edge magnetism. 
We point out that magnetic edge states might not exist 
in real systems. And even if they do, the intrinsic mag- 
netism would be so weak, that realistic applications in 
spintronic devices will not be feasible. In order to use a 
physical effect in a practical device, the effect must be 



(i) a strong intrinsic effect and (ii) it must be stable at 
room temperature. We show that graphene edge mag- 
netism does not fulfill these two criteria: (i) the effect is 
easily destroyed by a multitude of mechanisms and (ii) it 
is not stable at room temperature, even in perfect sam- 
ples. 



II. COMPUTATIONAL METHODS 

Our electronic structure calculations are based on den- 
sity functional theory (DFT) 17 , using a plane wave basis 
set (cutoff energy: 450 eV) and the projector augmented- 
wave (PAW) method 18 (we use the program VASP, ver- 
sion 4.6) j 19 i 20 The exchange-correlation interactions were 
approximated by the generalized gradient approxima- 
tion (GGA) within the Perdew-Wang parametrization 
(PW91)^i Total energies in the self-consistency cycle 
were converged such that energetic changes were less than 
10 -4 eV. The /c-space integrations were carried out us- 
ing the method of Methfessel and Paxton^ in first order, 
with a smearing width of 0.1 eV for the structural opti- 
mizations, and with the improved tetrahedron method 23 
for a final static calculation of the total energies. The 
optimal sizes of the k-point meshes for different systems 
were individually converged, such that changes in the to- 
tal energy were reduced to at least 3 meV/atom^ 

In our study we consider 10-ZGNRs and 12-ZGNRs, 
where the prefix '10' or '12' represents the number of 
zigzag lines across the width of the ribbon. A suffix '+H' 
indicates that each carbon edge atom is passivated by 
one hydrogen atom, '+EV that one edge contains an 
edge vacancy, L +zz(57y that the edge is reconstructed 
with alternating lines of pentagons and heptagons, and 
'+2211 ' that one out of three edge atoms is passivated 
by two hydrogen atoms and the other two by one hy- 
drogen atom. The atomic structures of some of these 
systems are shown in Figs. EEl For the calculations of 
a H2 molecule, graphene and the ZGNRs, we simulate 
free-standing objects, where the vacuum region between 
replica in neighboring unit cells is at least 10 A wide^ 

In Table H] we list the edge energies of all ZGNRs that 
are considered in this study. The edge energy quantifies 
the energy needed to form an edge from infinite graphene 
and (for hydrogen passivated systems) molecular hydro- 
gen, i.e. the enthalpy of the virtual reaction 

graphene + A H /2 H 2 — ► ZGNR. 

The edge energy can be defined per edge atom (E*^ ) or 
per unit length along the edge (E^ ) as follows 
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where Ac, Ah, and A c g are the number of carbon, hy- 
drogen, and carbon edge atoms per unit cell, respectively. 
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771 at 

-^edge 


ZTilen 
-^edge 


^ads 




(eV/at) 


(eV/A) 


(eV/at) 


10-ZGNR 


2.99 


1.21 


- 


10-ZGNR+H 


0.26 


0.11 


-2.73 


12-ZGNR 


3.00 


1.21 


- 


12-ZGNR+H 


0.28 


0.11 


-2.72 


12-ZGNR+Z2H 


0.08 


0.03 


-2.91 


12-ZGNR+^(57) 


2.39 


0.97 


- 


12-ZGNR+z*(57)+H 


0.87 


0.36 


-1.52 


10-ZGNR+EV 


2.88 


1.17 


- 


10-ZGNR+EV+H 


0.86 


0.35 


-2.02 



TABLE I. Edge energies E^ ge and i^dge as defined by 
equations [T] and O as well as hydrogen adsorption ener- 
gies E^ ds as denned by equation 03 both calculated at the 
DFT/GGA/PAW level of theory. 



the (magnetic) ground state as given in Tab. HH ^coh P 
= -9.244 eV/atom and E^ h = -3.394 eV/atom are the 
cohesive energies of a carbon atom in graphene and a hy- 
drogen atom in the H2 molecule, respectively. A is the 
lattice constant in the periodic direction, as indicated in 
Fig. Efb). Note that in Eq. [2] we divided by 2 A, because 
the unit cell contains both edges of the ZGNR. Our re- 
sults for edge energies compare well with previous data 
for such systems j 26 i 27 

The hydrogen adsorption energies in Tab. [I] quantify 
the gain in energy of a hydrogen passivated ZGNR, as 
compared to an unpassivated system and molecular hy- 
drogen. In other words, it is the enthalpy of the reaction 

ZGNR + A H /2 H 2 — > 'ZGNR+H'. 

The hydrogen adsorption energy per carbon edge atom 
is defined as 
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where £ t z ™ R+H and £ t z ? NR are the total energies of 
the hydrogen passivated and the non-passivated ZGNR, 
respectively. 

For magnetic systems we performed collinear, spin- 
polarized DFT calculations. The results are listed in 
Tab. HI comprising energies and magnetic moments for 
nonmagnetic (NM), ferromagnetic (FM), and antiferro- 
magnetic (AFM) states. The energies of different mag- 
netic states are compared via the quantity A£? mag , that 
is the difference in total energy with respect to the mag- 
netic ground state (E^) per magnetic atom: 
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is the total energy per unit cell of the ZGNRs in 



where Ama is the number of magnetic atoms per unit 
cell. For ZGNRs Ama is the number of magnetic carbon 
edge atoms per unit cell. 
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FIG. 2. (Color online) Atomic structures (top) and nonmag- 
netic band structures (bottom) of a 12-ZGNR with (a) unpas- 
sivated edges and (b) monohydrogenated edges. The fatness 
of the bands is proportional to the orbital character of the cor- 
responding wave function on the encircled edge atom (top). 
Orange (light gray) represents the n character of the bands, 
and blue (dark gray) represents their a character. Gray hor- 
izontal lines were added to the atomic structure (top), in or- 
der to indicate the boundaries of the corresponding unit cell 
along the periodic direction (see also Fig. [5]). Because of the 
presence of (partially) flat bands near Ef in the nonmagnetic 
band structure, there will be a magnetic transition into an 
antiferromagnetic ground state with a band gap (not shown) . 



III. DISCUSSION 

A. Stability of edge states 

In a planar sp 2 system like graphene, the electronic 
states split into in-plane (a) and out-of-plane (tt) states 
that are decoupled by symmetry. At an unpassivated 
zigzag edge the hexagonal carbon network is interrupted, 
and both the a and the tt system form edge states. The 
edge states of the a system are unpaired electrons in sp 2 
orbitals, i.e. dangling a bonds. The edge states of the 
7r system are the ones that were introduced above, and 
those are usually discussed in the literature. Within the 
nonmagnetic, electronic band structure of an unpassi- 
vated ZGNR both sets of edge states appear as nearly 
flat bands (or as flat parts of bands) at Ep 26 ^ 28 These 
bands are highlighted in Figs. [2]|2J In all of these plots, 
the fatness of the bands is proportional to the orbital 
character of the corresponding wave function at an edge 
atom. Orange (light gray) represents their tt charac- 
ter, and blue (dark gray) represents their a character. 



FIG. 3. (Color online) Atomic structures and band structures 
of a 12-ZGNR with a zz(57) edge reconstruction for (a) un- 
passivated edges, (b) monohydrogenated edges (see Fig. [2}. 
The system is metallic, nonmagnetic, and has edge states. 



The nonmagnetic band structure of an unpassivated 12- 
ZGNR is shown in Fig. Efa). Here the a and tt edge 
states exist in a range of about 0.5 eV above and below 
Ey- Only the degenerate part of the tt band has edge 
state character, and its non-degenerate part has bulk 
state character. Flat bands lead to a high DOS. As the 
chemical reactivity is proportional to the DOS near Ep, 
an unpassivated ZGNR would be highly reactive, since 
the edge states are available to form chemical bonds. In 
other words, the edge states do not only lead to a mag- 
netic instability, but they also lead to a chemical and 
structural instability. The latter will be discussed below. 
For the 12-ZGNR the magnetic instability leads to an 
AFM ground state with a band gap of 0.7 eV (Fig. [2ja) 
shows the band structure of the metastable NM state). 
The high instability of an unpassivated edge is marked 
by a rather high edge energy E^ of the 10- ZGNR and 

12-ZGNR systems of 1.21 eV/A (see Tab. IJ). 



1. Edge passivation 

One stabilization mechanism for an unpassivated 
ZGNR is edge passivation. The system will catch any 
bond partner it can possibly get to saturate its dangling 
bonds. The majority of papers on ZGNRs consider mono- 
hydrogenated edges with a single hydrogen atom per edge 
atom. In these systems the dangling a bonds are satu- 
rated, and the a edge states are removed from the vicinity 
of Ep [see Fig. |2fb)]. However, the tt edge states per- 
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FIG. 4. (Color online) Atomic structure and band structure 
of a 12-ZGNR with a 2211 edge (see Fig. [2}. The system does 
not have edge states, it is nonmagnetic, and a semiconductor 
with a band gap of 0.61 eV. 



sist and render the system magnetic, which results in an 
AFM ground state discussed above (the band structure 
of the AFM state is not shown). For the 10-ZGNR and 
12-ZGNR the edge energy drops from 1.21 eV/A to 0.11 
eV/A after passivation, which stabilizes these systems 
rather significantly. 

However, this type of edge passivation does not cor- 
respond to the thermodynamic ground state, as shown 
by Wassmann et at. on the basis of DFT/GGA 
calculations^! According to these authors the mono- 
hydrogenated zigzag edge (labeled Z\ in their paper), as 
well as the zz(57) edges discussed below, are only stable 
at extremely low hydrogen concentrations (P ^ 10 -20 
bar at room temperature). At realistic hydrogen pres- 
sures, the only stable zigzag edge should be the Z211 state 
shown in Fig. [4j with two mono- and one dihydrogenated 
sites. Our calculation on a 12-ZGNR with the ^11 edge 
confirm these findings, as the corresponding edge energy 
is only 0.03 eV/A, i.e. more than a factor of three lower 
than the monohydrogenated edges. The resulting system 
does not have any edge states. It is a semiconductor with 
a band gap of 0.61 eV, and it is nonmagnetic (see Fig.|4]). 
The Z211 edge is by far the most stable zigzag edge that 
is known up to now. 

For the ^11 passivation type the dihydrogenated edge 
atoms are sp 3 hybridized. Therefore a departure from the 
regular sp 2 structure at the edge will be a perfect way to 
get rid of the edge states, thus removing the peak within 



the DOS at E? and stabilizing the overall structure. 



2. Edge reconstruction 

Another stabilization mechanism for an unpassivated 
ZGNR is a structural deformation via edge reconstruc- 
tion. In their DFT/GGA calculations Koskinen et 
at. show that a planar reconstruction with alternating 
pentagon-heptagon carbon rings along the edge, named 
zz(57) (see Fig. [3j), is 0.35 eV/A lower in edge energy 
than the normal zigzag edge. 26 In our own calculations 
of a 12-ZGNR+zz(57) system we found that the recon- 
structed edge is 0.24 eV/A lower in edge energy than the 
12-ZGNR system (see Tab.[Ij). Compared to the edge en- 
ergies of the hydrogen passivated systems this is only a 
modest lowering, indicating that the zz(57) edge recon- 
struction is likely to occur only in high vacuum. 

The Koskinen reconstruction leads to the formation of 
triple bonds in the nearly linear parts of the heptagons, 
which shift the a edge states away from Ep [see Fig.[3fa)]. 
The 7r edge states remain near Ey in two normally dis- 
persed bands, that render the system metallic. However, 
there is no magnetic transition, because these bands are 
not flat. We confirmed this fact on the basis of spin- 
polarized calculations, which all converged into a NM 
state. Again, only the degenerate part of the band at E? 
can be associated with edge states, the non-degenerate 
part has bulk character. In the hydrogen passivated 12- 
ZGNR+zz(57)+H the a edge states are removed from 
the band structure, and the occupied tt bands are prac- 
tically unchanged [see Fig. [3^b)]. 

Because of the absence of a states near Ep Koskinen et 
al. argue that zz(b7) reconstructed edges do not need to 
be stabilized by passivation. They further show that such 
systems do not favor hydrogen passivation, and therefore 
they call the zz(57) reconstruction "self-passivating" . 
Our results support this judgment only partially. First, 
the edge energy of 0.97 eV/A for the unpassivated 12- 
ZGNR+£z(57) is still very high. And second, we find 
that the hydrogen adsorption energy is —1.52 eV per edge 
atom. Here a negative sign indicates that hydrogenation 
leads to an increase of binding energy, i.e. that hydrogen 
adsorption is favored. After passivation the edge energy 
of the system (12-ZGNR+zz(57)+H) drops to a value 
of 0.36 eV/A, indicating that the passivated system is 
energetically much more favorable than the unpassivated 
system. However, because of the absence of a states near 
Ep, hydrogenation is the least favorable process for these 
particular systems, in contrast to all other systems. This 
is reflected in a lower hydrogenation energy, and a rela- 
tively high edge energy of the 12-ZGNR+zz(57)+H sys- 
tem. 

Among other unusual edge reconstructions, the exis- 
tence of the zz(57) type was confirmed by analyzing 
aberration-corrected transmission electron microscopy 
(TEM) graphs of free-standing graphene i 29 i 3 ° i31 



System 


State 


iV at 


Ob) 


^tot 

(eV) 


(meV/at) 


10-ZGNR 


NM 

FM 

AFM 


60 
60 
60 


0.0 

7.7 
0.0 


-535.04379 
-536.62974 
-536.64359 


269 
4 



12-ZGNR 


NM 

FM 

AFM 


72 
72 
72 


0.0 

7.7 
0.0 


-645.97119 
-647.54796 
-647.56282 


265 
2 



10-ZGNR+H 


NM 

FM 

AFM 


66 
66 
66 


0.0 
1.4 
0.0 


-573.24160 
-573.36746 
-573.40631 


27 
6 



12-ZGNR+H 


NM 

FM 

AFM 


78 
78 
78 


0.0 
1.4 
0.0 


-684.05304 
-684.20262 
-684.22471 


29 
4 



10-ZGNR+EV 


NM 
FM 


59 
59 


0.0 
4.0 


-527.32015 
-528.09920 


260 
a 


10-ZGNR+EV+H 


NM 
FM 


64 
64 


0.0 
1.0 


-557.08600 
-557.18430 


33 

Q a 


Fe 


NM 
FM 


1 
1 


0.0 
2.1 


-7.7566476 
-8.1517557 


395 



NiO 


NM 

FM 

AFM 


12 
12 
12 


0.0 

4.8 
0.0 


-68.678312 
-68.722936 
-70.142254 


244 

237 





a Only the defect-free edge carries magnetic moments. 

TABLE II. Magnetic states of different systems calculated at 
the DFT/GGA/PAW level of theory. "State" denotes states 
with different magnetic ordering: NM (non-magnetic), FM 
(ferromagnetic), AFM (antiferromagnetic). N a t is the number 
of atoms per unit cell, fi is the total magnetic moment per 
unit cell, Et t is the total energy per unit cell. AE mag is 
denned in Eq. UJ 



3. Edge closure 

Finally, edge closure in mult i- layered graphene is the 
third stabilization mechanism we would like to empha- 
size. Liu et at. study the edge structures of graphite 
by atomically resolved high-resolution TEM. 32 In a se- 
ries of tilting experiments the authors demonstrate, that 
after thermal treatment (at 2000°C) the zigzag and arm- 
chair edges between adjacent graphene layers are mostly 
closed, and that they tend to form folded layers. Edge 
closure in multi-layered graphene therefore seems to be 
a simple way to get rid of edges, as well as of the corre- 
sponding edges states [see Fig. DJc)]. 



B. Stability of edge magnetism 

To summarize the first part of our discussion, we 
showed that ZGNRs with magnetic edge states are not 
very likely to exist. Nevertheless let us now hypothesize 
that such systems could be made, and let us discuss the 
resulting edge magnetism in more detail throughout the 
second part of this paper. 

We start our discussion with two standard magnetic 
systems, the FM iron (bcc phase, space group Fm3m) 
and the AFM NiO (rhombohedral lattice system, space 



group R3m). The energies and magnetic moments of 
different magnetic states for the two systems are given 
in Tab. HI For both systems the ground state mag- 
netic configuration is separated from other possible mag- 
netic configurations by an energy gap of several hundred 
meV per magnetic atom. This large energy gap is the 
reason, why magnetic order persists at room tempera- 
ture. The size of this energy gap AE^^ 1 , as given by 
the value of AE mag for the second lowest (GS+1) en- 
ergy state, is an indication of the thermal stability of the 
magnetic state. Ai^f^ 1 defines an an upper bound for 
the critical temperature T c max for spontaneous magneti- 
zation (Ai^f^ 1 = kT c max ). The real critical tempera- 
tures are usually much smaller. For iron and NiO the 
we find T c max to be about five times bigger than the ex- 
perimentally measured T c . Below we will use T c max to 
estimate the thermal stability of spontaneous magneti- 
zation in ZGNRs. 



1. Ideal ZGNRs 

Now consider the magnetic states of the ideal, unpas- 
sivated and monohydrogenated systems 10-ZGNR, 12- 
ZGNR 10-ZGNR+H, 12-ZGNR+H. The 12 carbon chain 
systems are shown in Fig. O Table [XT] indicates that for 
each of these systems a NM, FM, and AFM state can be 
found. For the FM configuration all magnetic atoms are 
aligned along the magnetic axis, and in the AFM con- 
figuration the atoms are ferromagnetically ordered along 
one edge, and antiferromagnetically ordered between op- 
posite edges. Figure [5fa) depicts the spin density of the 
AFM state of the 12-ZGNR. All of these results are in 
excellent agreement with other DFT studies of similar 
systems. 33 35 

The AFM configuration is the lowest energy state of 
all systems. However, the FM states are only AE^^ 1 = 
2 — 6 meV per edge atom higher in energy. The up- 
per bound for the critical temperatures in these cases is 
jimax = 70 K. Furthermore, these very small energy dif- 
ferences are of the order of the precision of modern DFT 
calculations. So within DFT the AFM and FM states can 
hardly be distinguished at all from an energetic point of 
view. The reason for this finding is the large width of the 
ribbons (19.8 and 24.1 A for 10-ZGNR and 12-ZGNR, re- 
spectively), which leads to a decoupling of the magnetic 
edges. This is in perfect agreement with the well known 
finding, that for large ribbon widths, the FM and AFM 
states are degenerated For ZGNRs with smaller widths 
the interactions between the tails of the edge states are 
more pronounced, and the energy differences between the 
AFM and FM configurations Ai^j^+i increase. It has 
been shown by Jung et al that A£^f + 1 follows a C/W 2 
law, where W is the ribbon width and C = 2.7 eVA 2 . 36 
But even for very narrow ZGNRs, AE^ + 1 will not sig- 
nificantly exceed 30 meV. 

We conclude that the edge magnetism in ideal, unpas- 
sivated and monohydrogenated ZGNRs will not be stable 



(a)12-ZGNR 




(b)10-ZGNR+EV 




FIG. 5. (Color online) The spin density of (a) an unpassi- 
vated 12-ZGNR in the AFM state and (b) an unpassivated 
10-ZGNR with an edge vacancy in the FM state. The red 
(dark gray) contours on the right hand side represent posi- 
tive spin densities within the graphene layer; the yellow (light 
gray) contours on the left hand side represent negative den- 
sities. A is the lattice constant along the periodic direction. 
(b) Due to the presence of the defect, there are no remaining 
magnetic moments on the left edge. 



at room temperature. Therefore practical applications in 
spintronic devices will not be feasible. We estimate that 
the AFM state can only be stabilized at temperatures 
T < 10 K. Above that temperature ideal ZGNRs will 
simply exhibit paramagnetic behavior. 

These findings are in good agreement with a recent ex- 
perimental study by Sepioni et al£L Those authors stud- 
ied graphene laminates by SQUID magnetometry, and 
found that the laminates are strongly diamagnetic and 
exhibit no sign of ferromagnetism down to 2 K. Below 
50 K they detect a very weak paramagnetic contribution. 
These results can be explained by our previous argument, 
that the majority of edges is either reconstructed, passi- 
vated, or closed. Only a few remaining magnetic edges 
would give rise to the very weak paramagnetic contri- 
bution. This explanation, however, cannot account for 
very narrow distribution of magnetic moments that was 
measured in Ref. [37j. 

A further confirmation of the highly unstable mag- 
netism in ZGNRs was given by Yazyev et al. 38 They used 
a DFT parametrized spin-Heisenberg model to show, 
that the spin correlation length of the magnetic edges 
is only about 1 nm at room temperature. 



2. Edge Defects 

Under current experimental conditions one cannot pro- 
duce such ideal ZGNR. It is also rather questionable, 
whether this will be possible at all. In reality one always 
has to deal with non-ideal structures, containing impuri- 



ties, vacancies, lattice defects, etc. 39 

Figure [5fb) shows a 10-ZGNR with an edge vacancy 
(EV) on the left hand side of the ribbon. This vacancy 
totally suppresses any magnetic moments on that side, 
but it does not influence the moments on the other side. 
Therefore only one magnetic solution (FM) exists for 10- 
ZGNR+EV and 10-ZGNR+EV+H. Apart from that, the 
values of AE mag for the FM an NM states are similar to 
previous cases without edge vacancy. This result shows 
again that opposite edges are almost decoupled. 

For the unpassivated 10-ZGNR+EV A£gf + 1 is very 
high. This is an effect of the dangling a bonds, which 
have unpaired electrons with individual magnetic mo- 
ments of about 1 /iB • These unpaired electrons render the 
system strongly magnetic. However, as discussed above, 
unpassivated systems are highly unstable and they will 
never exist at ambient conditions. The instability of this 
system is again reflected by its high edge energy of 1.17 
eV/A (see Tab. HJ. It is interesting to notice that the 
edge energy of the 10-ZGNR+EV is lower than the edge 
energy of an ideal 10-ZGNR by 0.04 eV/A. Thus for un- 
passivated systems the formation of EV is energetically 
favorable. The reason is that the EV removes a edge 
states from the vicinity of Ep- This finding might be 
important for the understanding of disorder during the 
growth process, where unpassivated edges might exist at 
an early stage. 

However, if the edges are passivated, as for the 10- 
ZGNR+EV+H system, no a edge states are present. 
Consequently the formation of EV will not be favored any 

more (E^ ge is °- 21 eV M ni g her than for 10-ZGNR+H). 
As the 10-ZGNR+EV+H system has magnetic moments 
only along one edge AE^&g 1 = 33 meV can be consid- 
ered as the energy of moment formation along a single 
edge. The corresponding upper bound for the critical 
temperature T c max = 380 K is significantly higher than 
for ideal ZGNRs (T c max = 70 K). However we want to re- 
mind the reader that real critical temperatures are usu- 
ally much lower than T c max . So room temperature appli- 
cations of edge magnetism along a single edge are also 
unlikely. 

Huang et al. systematically studied the effect of edge 
defects (vacancies and substitutional dopants) on the 
magnetism with DFT/LDA calculations^ The stability 
of the magnetism is found to continuously decrease with 
increasing concentration of the defects. ZGNRs become 
nonmagnetic at concentrations of about one edge defect 
impurity per 10 A (in our example the concentration is 
one defect per 7.4 A). The authors judge that such critical 
defect concentrations are accessible in real samples. The 
corresponding suppression of edge magnetism is mainly 
caused by the removal of edge states from E^. 

Yazyev et al. used DFT/GGA calculations to show 
that the magnetic domain wall creation energy of 114 
meV for the ideal zigzag edge decreases dramatically in 
the presence of edges defects, down to energies of the 
order of 30 meV. This means that magnetic alignment 
along one edge, which might be relatively stable in ideal 




FIG. 6. The effect of charge doping in 12-ZGNR+H: AN 
is the number of excess electrons per edge atom, \i is the 
magnetic moment per edge atom in the FM state, and AEp 
is the change of the Fermi energy with respect to the undoped 
case for both the NM and the FM state. 



ZGNRs, will easily be destroyed by the occurrence of edge 
defects^ 

We therefore conclude that the edge magnetism, which 
is already a very weak effect in ideal ZGNRs, is further 
weakened or even totally disappears in the presence of 
edge defects. The stability of the edge magnetism could 
only be enhanced, if high density edge defects are con- 
fined to one edge of the ribbon, while the other edge 
remains in an ideal zigzag shape (10-ZGNR+EV+H). 
However, even in this very artificial case, the magnetism 
is still too weak to allow for room temperature applica- 
tions. 



the band structure of the undoped 12-ZGNR+H system 
[see Fig. 12(b)] is symmetric near .Ep, and it is becoming 
more and more asymmetric away from Ep. Furthermore 
for large doping levels the rigid band model breaks down, 
and the band structure of the doped system changes sig- 
nificantly, as compared to the undoped case. 

The critical doping level for a 12-ZGNR-H is ca. 0.5 
electrons per edge atom, which agrees very well with 
previous results j 41 i 42 However this is a very high doping 
level, which is not reached in standard experimental set- 
tings. The charge fluctuations in typical graphene charge 
puddles are of the order of n^D = 10 11 cm -2 , and back- 
gate doping can inject charge carriers in concentrations 
up to n 2 D = 10 12 - 10 13 cm" 2 ^ 3 -! 4 ^ Such substrate in- 



duced charge concentrations {n^D 



10 



li 



10 



13 



cm 



l ) 



would correspond to doping levels AN between 0.0006 
and 0.06 electrons/holes per edge atom for the 12-ZGNR- 
H ribbon. 45 Under such conditions the weakening of the 
edge magnetism would only be marginal. 

Figure [6] allows for the translation of the critical doping 
level into a critical shift in Ep, which is about ±2 eV. So 
by contacting the 12-ZGNR-H at two points along the 
periodic direction, and by applying critical voltage of 4 
eV (shifting the electronic structure by +2 eV on one 
contact and by -2 eV on the other contact), one would 
also destroy the edge magnetism close to the contacts. 
Smaller voltages may not destroy the edge magnetism, 
but it will be weakened quite significantly. As the charge 
mobility in graphene is very high, voltages of the order 
of a few volts can indeed be applied 46 , which leads to a 
significant weakening of the edge magnetism. 

We thus conclude that under current experimental con- 
ditions, back-gate doping and charge puddles do not sig- 
nificantly influence the edge magnetism in ideal ZGNRs 
of rather small width. However, when contacting ZGNR 
electrically, and after applying voltages of a few volts, the 
edge magnetism will be significantly weakened. 



3. Charge 



In experimental measurements of graphene systems 
one always has to deal with (intentional or uninten- 
tional) charge doping effects, which are induced by the 
substrate/back-gate or by various impurities. 

In Fig. [6] we show the effect of charge doping on a 12- 
ZGNR-H system. AiV is the number of excess electrons 
per edge atom. Positive values represent electron doping, 
and negative values represent hole doping. We normalize 
A TV to the number of edge atoms, because the edge states 
are the states, which are directly affected by the doping. 
We see that progressive charge doping by electrons or 
holes move Ep away from the aforementioned peak in 
the DOS (at AE-p = 0), which gradually suppresses the 
magnetic transition. Therefore the magnetic moment per 
edge atom (/i FM ) becomes smaller and smaller until it 
vanishes eventually, and the corresponding system be- 
comes NM. For small values of ±A7V electron and hole 
doping is symmetric, but for larger values the two scenar- 
ios become more and more asymmetric. This is because 



IV. CONCLUSIONS 

In summary, we critically discussed the stability of 
edge states and edge magnetism in ZGNRs. In the 
first part of our paper we pointed out that magnetic 
edge states are not very likely to exist. We showed 
that there are at least three very natural mechanisms 

- edge reconstruction, edge passivation, and edge closure 

- which dramatically reduce the effect of edge states in 
ZGNRs, or even eliminate them completely. In the sec- 
ond part of our paper, we showed that if ZGNRs with 
magnetic edge states could be made, the intrinsic mag- 
netism would not be stable at room temperature. Fur- 
thermore, charge doping and the presence of edge de- 
fects further destabilize this already rather weak type of 
collective magnetism. We conclude that antiferromag- 
netic ZGNRs might only be observed under ultra-clean, 
low-temperature conditions in defect-free samples. The 
present state of research indicates that edge magnetism 



within graphene ZGNRs is much too weak to be of practi- 
cal significance, in particular for spintronics applications. 
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